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We discuss the problem of the bias of the Internal Linear Combination (ILC) CMB map and 
show that it is closely related to the coefHcient of cross-correlation K{£) of the true CMB and 
the foreground for each multipole £. We present analysis of the cross-correlation for the WMAP 
ILC quadrupole and octupole from the first (ILC (I)) and the third (ILC (III)) year data releases 
and show that these correlations are ^ —0.52 — 0.6. Analysing lO** Monte Carlo simulations of 
the random Gaussian CMB signals, we show that the distribution function for the corresponding 
coefficient of the cross-correlation has a polynomial shape P{K,£) oc (1 — K'^Y~^. We show that 
the most probable value of the cross-correlation coefficient of the ILC and foreground quadrupole 
has two extrema at iiT ~ ±0.58. Thus, the ILC(III) quadrupole represents the most probable 
value of the coefficient K. We analyze the problem of debiasing of the ILC CMB and pointed out 
that reconstruction of the bias seems to be very problematic due to statistical uncertainties. In 
addition, instability of the debiasing illuminates itself for the quadrupole and octupole components 
through the flip-effect, when the even £ + m modes can be reconstructed with significant error. 
This error manifests itself as opposite, in respect to the true sign of even low multipole modes, 
and leads to significant changes of the coefficient of cross-correlation with the foreground. We 
show that the CMB realizations, whose the sign of quadrupole (2, 0) component is negative (and 
the same, as for all the foregrounds), the corresponding probability to get the positive sign after 
implementation of the ILC method is about 40%. 
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1. INTRODUCTION 

Since the COBE experiment and then af- 
ter the Wilkinson Microwave Anisotropy Probe 
(WMAP) [D null [3 El [7] first and third year 
data releases, the problem of reconstruction of the 
CMB low multipoles < 10 (including the power 
of a quadrupole, planarity and alignment of the 
quadrupole and the octupole etc.) attracts very seri- 
ous attention. A lot of cosmological models have been 
involved for the explanation of the peculiarities of the 
ILC (Internal Linear Combination ) low multipoles 
including the model of running spectral index of the 
primordial adiabatic perturbations [11 El [3l H [5] , non- 
trivial topology of the Universe [SI [S] , Broken Scale 
Invariance (BSI) of the power spectrum [TUl [H] , pri- 
mordial magnetic field [HI [131 111] i Bianchi Vllh cos- 
mological model [m Hg nil nil US] , etc. 

The problem of low power of the CMB quadrupole 



can be related to the non-Gaussianity of the WMAP 
low multipoles, including planarity and alignment of 
the quadrupole and octupole [2Qj ^2T, '22', "23] , cross- 
correlations with the foregrounds ^24, 25j 2^1 HI [?7]. 
local defects of the CMB map [28] , etc. These pecu- 
liarities need some theoretical explanation, since the 
origin of those non-Gaussian features still remains un- 
known. 

Another way to look at the problem of the low 
multipole range of the ILC CMB signal is to take 
under consideration the cross-correlation between the 
CMB signal and the foregrounds, known as a bias 
of the ILC map [29l|4l [27l [24l ISQ], some effects of 
systematics [31] including contamination of the £ = 2 
mode by dipole through far sidelobes [35]. 

In this paper we would like to re-examine one of 
the most important issues of the data analysis of the 
modern CMB experiments related to the problem of 
the CMB~foreground separation by the ILC Method. 
Needless to say, that separation of the CMB signal 
and foregrounds is the major problem for all the CMB 
experiments including the ongoing WMAP and the 
future Planck missions. The basic formalizm of the 
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ILC method which aUows one to reconstruct the CMB 
signal from the observational data sets is to combine 
the data in one (Linear Combination) map by using 
real weighting coefficients. Then, to find these weight- 
ing coefficients, we need to minimize some functional 
defined on multifrequency CMB data sets containing 
information about the primordial CMB signal, the 
instrumental noise and the foregrounds. This func- 
tional (particularly, variance by Tegmark and Efs- 
tathiou [33j (hereafter TE)) can be defined both in 
the pixels domain (like it was done by the WMAP 
team), and in the spherical harmonics domain [34] 
(hereafter TOH), and [1 [30l [35] . Below, we will dis- 
cuss only the low CMB multipoles {£ < 10) and, par- 
ticularly, the quadrupole for which the ILC approach 
in multipole domain seems to be very useful. 

However, the ILC approach in combination with 
the Minimal Variance Method (hereafter, ILC MVM) 
produces a specific feature of the derived CMB sig- 
nal, known as the bias. The bias reflects directly the 
cross-correlation with foregrounds, and leads to the 
non-Gaussianity and statistical anisotropy of the de- 
rived CMB through coupling with the foreground. 
This is because the MVM was originally designed 
for separation of the Gaussian signals, of which the 
statistical properties are completely specified by the 
variance of each component. At the same time, it is 
well known in the data analysis that the MVM is 
not motivated properly, especially in cases when we 
have only a unique realization of the random process, 
where the CMB signal is mixed with the significantly 
non-Gaussian foreground. Moreover, the statistics of 
the low multipoles of the CMB signal are extremely 
peculiar, since we do not have the statistical ensemble 
of realizations. 

In general, this paper is devoted to investigation 
of the ILC bias and its statistical uncertainties. We 
will show that so-called debiasing of the ILC CMB 
has potential instabilities and unpredictability for a 
single realization of the CMB. 

The outline of the paper is as follows. In Section 2, 
we focus our attention on mathematical basis of the 
ILC MVM exactly in the same way as it was proposed 
by [1]. By analysis of the model with homogeneous 
foreground spectra, we show that the bias of the ILC 
image, and the bias of the power spectrum are in or- 
der of K{i) and K^{£), where K{£) is the coefficient of 
"the true CMB — true foreground" cross-correlation. 
Following [26] , we will call this cross-correlation "the 
cosmic covariance". Because of statistical origin of 
this correlation, the cosmic covariance is uncertain 
for each single realization of the CMB signal, even 
for well known foregrounds. At the same time, de- 
rived by the ILC MVM, the ILC CMB signal and cor- 
responding ILC MVM foregrounds have zero cross- 



correlations due to design of the ILC MVM. 

In Section 3, we draw our attention to the investi- 
gation of the statistical properties of the coefficient of 
cross-correlation between the random Gaussian CMB 
and the WMAP foregrounds for K-W bands. We gen- 
erate 10^ realizations of the Gaussian CMB signals for 
the WMAP best fit ACDM cosmological model, and 
cross-correlate them with the foregrounds obtained 
by the WMAP team (the Maximum Entropy Method 
(MEM) foregrounds, which constitute as a sum of 
the synchrotron, free-free and dust emission for each 
frequency band). We will show, that the probabil- 
ity distribution function for these cross-correlations 
has a polynomial shape P{K,£) oc (1 — K'^Y~^. Af- 
ter that, we will find the most probable values of 
the coefficients K{£) and K {£). In another words, 
these coefficients determine the most probable value 
of the bias for the ILC CMB and the power spectrum. 
Taking under consideration that the maps ILC (III) 
of three years and ILC(V) of five years data are 
already debiased, as the WMAP team claimed, we 
check out the cross-correlation between the ILC and 
the WMAP forground. After debiasing, these cross- 
correlations should represent the possible value of the 
K{£) for the true CMB and the true foreground. For 
the WMAP K-W bands, the corresponding coeffi- 
cients of the cross-correlations are about K ~ —0.51. 
Remarkably, this value is close to the most probable 
value K{£) of the distribution function KP{K, £ — 2). 

In Section 4 we discuss some peculiarities of the 
ILC(III) quadrupole and the octupole. As it fol- 
lows from the Monte Carlo simulations presented 
in [3^ for low multipoles £ — 2,3, the sign of 
the {£ =2,™ =0), (2,2) and (3,1) components was 
changed to the opposite one for about 20% of the out- 
put CMB maps. This fiip-effect significantly changed 
the morphology of the output CMB map. We will 
show that the fiip-effect takes place for those input 
CMB quadrupoles, for which the sign of ^ = 2, to = 
mode is opposite to the sign of the foreground. For 
the Gaussian CMB signal about 50% of the realiza- 
tions have sign C2.0 ~ —sign i^2,o, where C2,o and i^2,o 
are the m=0 components of the quadrupole for the in- 
put CMB and the foreground, respectively. Roughly 
for 40% of those realizations, reconstruction of the 
^ = 2, TO = component could have a wrong sign of 
the (2,0) mode. In this Section we present some ana- 
lytical description of the flip-effect based on the prop- 
erties of the bias. 

We summarize our results in Conclusion. 
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2. MATHEMATICAL BASIS OF THE 
ILC MVM METHOD 

Let us draw our attention to the general formalism of 
the ILC method, assuming that we have experimental 
data AT{6p, 4>p,iyi) for the CMB anisotropy obtained 
for some frequency bands i^i , ■ • • ■ Here Bp and (f)p are 
the standard polar and azimuthal angles of the spher- 
ical coordinate system and the index p mark a cor- 
responding pixel. The observed sky temperature for 
each pixel p in a given band i is the linear combination 
of the true CMB signal Tc{p) and the true foreground 
fi (p) . In reality, the components of the foreground at 
different frequency bands (mainly the synchrotron, 
the free-free and the dust emission) have angular vari- 
ations of the spectral indicies. However, by using dif- 
ferent disjoint regions of the sky with nearly uniform 
foreground spectra and low level of spatial variations 
we can simplify the mathematical basis of the ILC 
method [4]. Following the WMAP model of the fore- 
grounds, let us assume that for some region of the 
sky the spatial variation of the foregrounds emission 
are negligible and fi{p) — SiF{p), where Si is the fre- 
quency spectrum and F{p) is the spatial distribution 
of the foreground for each band i. The main idea of 
the ILC method is to estimate the CMB signal by 
using real weighting coefficients Wi for each region of 
the sky with monotonic spectrum of the foreground 



Tiic{p) 



^T^{p) = T,{p)+VF{p),Y,W^ = U^) 



where F = WiSi. The coefficients Wi can be found 
by minimization of the variance (see [4]) 



Here is the varianc( 
variance of the true foreground and 



(2) 



Here is the variance of the true CMB, cr|, is the 



cTcF = {Tc{p)F{p)) - {T,{p)){F{p)) 

is the cosmic covariance between the true CMB and 
the foreground [26j . The angle brackets in Eql2] and 
hereafter denote the averages over the pixels belong- 
ing to each zone of the map with the monotonic spec- 
trum of the foreground. The minimum of the variance 
afi^ can be easily found from the equation -jS^ — 0, 
whence we get [J] 

Tiicip) = Tc{p) 



dwi 



^F{p), a 

F 



(J7 



2 

ilc 



(3) 



Let us define the coefficient of cross-correlation be- 
tween the true CMB signal and the true foreground as 



(4) 



Then, after substitution of Eq.lgl) in Eq.(I3]) we 
have 



Tuc{p)^T,{p)-K^F(p), al. 



OF 



^c(l-«')- (5) 



As one can see from Eq. ^ , the temperature per pixel 
Tiicip) the variance af^^ are biased. For Tudp) the 
difference Tiic{p) — Tc{p) is proportional to the k, and 
can be both positive and negative. For the variance 
(jfic, the bias is always negative and proportional to 
1 — . The existence of such a bias is not suppris- 
ing. The shift of the ILC variance and corresponding 
power spectrum was discussed in [4j [29l [27l [26] . How- 
ever, it would be important to note that the bias of 
the variance and the power spectrum is determined 
by the cosmic covariance, which is unknown for each 
single realization of the CMB sky. This is the reason 
why it can not be corrected by any methods because 
of the statistical uncertainties. 

The ILC approach presented above in the pixel 
domain can be easily generalized for the multi- 
pole domain by spherical harmonics decomposition: 

oo I 

e=o m=-l 

where \aim\ and (\>i^m are the moduli and phases of 
the coefficients of the expansion. Direct conversion of 
Eq.([5|) to the multipoles coefficients gives us 



-ilc 



-Fc 



OF 



(7) 



It is easy to show that this solution for the derived 
ILC signal c^'^ can be obtained by minimization of 
the variance 



3 |2 



(8) 



Moreover, follow to TOH approach we can define the 
ILC CMB in multipole domain by using weighting 
coefficients Wi{£),J2i^ii^) — 1 and minimizing the 
power spectrum for each multipole 



(i) |2 



For the model of uniform spectra of the foreground 
this approach leads to 



(El=-f |c£,mP)^ „(2) 
--fa 



\F^ ' 



(9) 
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Table 1: The table of Kff{£) coefficients for "V-W", "Ka-V", "Ka-W" and "Q~V" MEM foregrounds. The 



first column indicates the multipole £. 



Kffii) 




Ka-V 


Ka-W 


Q-V 


2 


0.999969 


0.999014 


0.998761 


0.999481 


3 


0.994699 


0.999534 


0.995187 


0.999770 


4 


0.996702 


0.999497 


0.996098 


0.999708 


5 


0.994615 


0.998862 


0.995825 


0.999647 


6 


0.995481 


0.999354 


0.993799 


0.999563 


7 


0.993245 


0.998620 


0.995323 


0.999506 


8 


0.994281 


0.999252 


0.993458 


0.999575 


9 


0.993532 


0.998781 


0.994617 


0.999571 


10 


0.994638 


0.999068 


0.993319 


0.999463 



where K{£) is the coefficient of cross-correlation be- 
tween the ILC CMB and the foreground: 



K{£)^ 



^ilc 12 
rn I i,rn I 



(10) 



Taking Eq.® and Eq. (fTO|l into account, one can get 



C'^^£)^Cci£) [1-K\£)] 



(11) 



Before any future analysis, let us draw our atten- 
tion to the applicability of the model of the uniform 
foreground spectra, taking the MEM foregrounds in 
K-W bands of the WMAP data into account. For 
these foregrounds (sum of the synchrotron, free-free 
and dust emission for each frequency band) the coeffi- 
cient of cross-correlation is defined in the same way as 
Eq.lUni) with substitution of F^/l, F^^^l for 4'^, 
inEq.Hni). 

In Table[T] we show the coefficients of cross- 
correlation for different combination of the WMAP 
MEM foregrounds. One can see that corresponding 
coefficients of cross-correlation between different fore- 
grounds (for different bands) are extremely close to 
unity. Thus the model of uniform spectra of the fore- 
ground is very well motivated. As it follows from 
Eq.® and Eq. fTTj) . the corresponding bias of the ILC 
coefficient is proportional to the K{£), while the 
bias of the power spectrum is of the order of K'^{£). 
Since the true CMB constitute as a random (Gaussian 
?) process, the coefficient of cross-correlation K{£) 
represents a random (non-Gaussian) process as well. 
Thus, for each particular realization of the CMB sky, 
the coefficient of cross-correlation K{£) remains un- 
certain. Note that the WMAP team, as mentioned 
in [57], to correct the CMB power obtained by the 
ILC method took under consideration 10^ 10'* real- 
izations of the random Gaussian CMB, creating the 
statistical ensemble of realizations. However, even af- 
ter the averaging over realizations we can not predict 



the bias for one particular realization of the CMB due 
to the cosmic covariance. Moreover, as it is seen from 
Eqs.® and (Hi]), the average of the ILC CMB sig- 
nal over the ensamble of realization reduces the bias 
factor (K) down to zero, while {K'^) ^ 0. Thus, for 
the ensemble of realization the correction of the CMB 
power can be done successfully, while the correction 
of the signal itself seems to be very problematic. 

One more problem is related to the statistical 
properties of the true CMB signal. As it is well known, 
the low multipoles {£ < 10) of the WMAP CMB sig- 
nal reveal significant departure from the statistical 
isotropy and homogenity 36, 37J. If this effect is the 
manifestation of the primordial non-Gaussianity of 
the CMB, we simply do not have the correct model of 
the statistical ensemble of realizations. This remark 
illuminates significantly stronger the problem of in- 
stability of the reconstruction of the low multipoles 
of the CMB signal. 

Not to be able to give the general solution of 
the problems mentioned above, we would like to pro- 
pose indirect method of evaluation of the bias in the 
image and power domains by looking at the cross- 
correlation coefficient K{£) of the derived CMB sig- 
nal and the MEM foregrounds, including the the syn- 
chrotron map of Haslam et al. [35] (hereafter HA). 
The idea is based on the fact that in Eq. (O the cross- 
correlation of the ILC and the foreground is exactly 
zero. Thus, if by implementation some methods of de- 
biasing of the ILC map we will get K ^ 0, it would 
be reasonable accurate estimator of the initial coeffi- 
cient of cosmic covariance k. In next section we will 
discuss this approach in details. 

3. THE ILC(I) AND ILC(III) CROSS- 
CORRELATIONS WITH THE 
FOREGROUNDS 

As was pointed out by the WMAP team [4l the cor- 
rection of the ILC CMB can be done by implementa- 
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Figure 1: The quadrupole component of the ILC(I) (top left), the ILC(III) quadrupole (middle left) and the 
difference between them (bottom left). Right column is for the ILC(I), ILC(III) and difference of octupole. 



tion of additional information about the foregrounds 
and 100 realizations of the random Gaussian CMB 
signals, which mimic the WMAP K-W bands in com- 
bination with MEM foreground. Repeating the ILC 
MVM method for each from these 100 realization, 
the WMAP team took under consideration some sys- 
tematic of the ILC method and made corresponding 
correction of the derived the ILC (III) map. 

It is important to note that the WMAP team [J] 
pointed out that difference between the ILC (III) and 
the ILC (I) is mainly due to the bias. In Fig.[T] we 
show the quadrupole and the octupole components 
of the ILC (I) and the ILC (III) maps and difference 
between them to illustrate the correction of the bias 
made by the WMAP team for low multipoles range 
of the signal. 

Let us draw our attention to the cross-correlation 
of the ILC (I) and ILC (III) maps with the MEM fore- 
grounds (being the sum over synchrotron, free-free 
and dust emission maps). In Fig. [5] we show the co- 
efficient of the ILC-foreground cross-correlation K{i) 
for the multipoles ^ = 2-10 and for the WMAP MEM 
foreground for K-W bancf^. In addition in Fig. [21 we 
show Kj{t) for the ILC(I), the ILC(III), and the 



^ For the ILC{I) we use the WMAP MEM foregrounds from 
the first year data release. 



Haslam et al. synchrotron map [38]. As one can see 
from this figure, the quadrupole {i = 2) and £ = 5 
multipole of the ILC (III) have negative correlation 
{K ^ —0.5) with Ka-W foregrounds. Moreover, the 
ILC(III) quadrupole has K ^ —0.6 with the Haslam 
et al. quadrupole. 

3.1. Statistical properties of the bias 

To understand the properties of the ILC (III) — 
foreground cross-correlation, we performed numerical 
test taking under consideration 10000 realizations of 
the random Gaussian CMB maps and cross-correlate 
them with the same model of the foregrounds as in 
Fig. [2] without any debiasing. 

In Fig. [3] we show the probability density function 
P{Kj{£)) versus Kj{i) for the multipoles ^ = 2 10. 

For the input random Gaussian CMB signal the 
shape of the distribution functions is perfectly fitted 
by the function 

P{K,e)^Aiil-Ky-\ (12) 

where is the normalization constant. Using 
P{K, £) for the input signal shown in Fig. [3] (top left), 
we have found the first < K >= —0.00043 and the 
second < K'^ >= 0.19934 moments of PiK,£ = 2), 
which is in agreement with corresponding value < 
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/ \ 
/ \ 



Figure 2: Top. The coefficient of cross-correlation of the ILC(I) and ILC(III) with the foregrounds. The dash 
line is for the ILC (HI) and the Haslam et al. map (1982); the triple dots dash line is for the ILC(I) and 
Haslam et al. map; the thick solid line with stars is for the cross- correlation of ILC(I) and V-band MEM 
foreground. The dash line with diamonds is for the ILC(III) and the V-band MEM foreground. Bottom. The 
coefficients of cross- correlations between V- and W-foregrounds of the WMAP (the thick solid line), V- and 
Haslam et al. 1982 (the dash line) and finally ILC(I) and ILC(III) (the triple dots-dash line). The dash line 
is for Q- and V-bands, the stars are for Ka-V and the boxes are for Q-V foregrounds. 



Table 2: The table of S coefficients for 10000 realizations of the CMB input and output maps. The top row 
contains data for the quadrupole component, the bottom one does for the octupole. and 3 mark real and 
imaginary parts of the a^_„i coefficients. 



1 = 2 


2,0 


5R(2,1) 


5(2,1) 


5R(2,2) 


9(2,2) 






S = -\ 


2148 


361 


55 


1756 


257 






t = Z 


(3,0) 


3fi(3,l) 


3(3,1) 


3fi(3,2) 


5(3,2) 


5R(3,3) 


3(3,3) 


S = -\ 


232 


1852 


337 


106 


367 


533 


355 



K'^ >= 0.2 from Eq.p^. 

At the end of this section we would like to present 
the distribution functions for the first and second mo- 
ments of the coefficient of cross-correlation between 
the Gaussian realizations of the CMB signal and the 
foreground. These moments determine the bias of the 
ILC signal and the power spectrum. It is not suppris- 
ing that the first, KP{K) and the second K'^P{K) 
moments shown in Fig.lH have two points of extrema, 



which can be easily estimated from Eq. fT^ : 

where n is the order of moment K"P{K), and i is 
the multipole number. For £ = 2 and n = 1 we 
have K^^^ = ±0.577 and K^'^ = ±0.707. This re- 
sult clearly demonstrates that if the sign of k co- 
efficient is fixed, the most probable values of K is 
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Figure 3: The distribution function P{K) for cross-correlation between random CMB and the V-band fore- 
ground (left). From the top left down to the bottom right I — 2, 3. .10. J^ote that the shape of P{K) is the same 
for any foreground from K-W band including Haslam et al. (1982) signal. 



K^'^, and the most probable vale of is K^'"^. The 
bottom plot of Fig.|4] shows the value of the mo- 
ments M{n) = {K'^P{K,£)) for the ensemble of real- 
ization of the random Gaussian process. For the sec- 
ond moment n = 2 the shape of the function M{n) is 
{2£ + in agreement with estimation [27 of the 

bias of the power spectrum. From Eq. (jl3[) one can 
find the most probable value of for the quadrupole 
component K^^^f^ — 0.5. This means that the most 
probable value of the factor (1 — i^p^ob) ^-^^ ^^'^ 
then from Eq.dTTl), we get cdi = 2) = 2Cucie = 2). 
At the same time the most probable value of the 
bias of the ILC CMB is determined by the param- 
eter K — K^'^ = ±0.577. One can see that this value 
is remarkably close to the WMAP ILC (III) coefficient 
of the cross-corelation with the foreground. 



4. PECULIARITIES OF 

QUADRUPOLE END OCTUPOLE 

In addition to the cross-correlation with the fore- 
ground, wc have discovered one more feature of the 
ILC method, which could have a significant im- 
pact on solution of the problem of the ILC (III) 
quadrupole [39]. We call this effect as a flip of the 
sign of the (2,0) component of the ILC quadrupole 
with respect to the sign of the true CMB (2,0) com- 
ponent. To understand quantitatively this effect, let 
us come back to Eqs.(j9l) and (fTO| for the quadrupole 
component £ = 2 and take under consideration the 
(2, TO = 0) mode: 

|c2,ol cos(6,o) = |c2,o|cos(7y2,o) 



|F2,o|cOs((^2,o), 



(14) 



where |..| denotes the modulus, and ^2,07^72,0 and 
^2,0 are the phases of 2, component of the ILC 
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Figure 4: The distribution functions KP{K) (top) and K^P{K) (middle) for cross-correlation between random 
realizations of the CMB signal and the V hand foreground. Bottom plot shows the distribution of {K"'{i)) for 
£ = 2 — 10. The solid line corresponds to {K'^{£)), the dash line is for {K*{£))'^ , the stars are for {K(£)), the 
rhombuses are for {K^{£)). 



quadrupole, the true quadrupole, and the foreground 
correspondingly. Since the m = mode has the 
only real part, these phases simply mean the sign 
of the (2,0) components. For our further analysis it 
is important to note that for all K~W foregrounds 
|-p2,o| > I-F24I, |-F2,2| and 02, = tt. Taking into ac- 
count Eq. (fTO|l . one can easily find that 



^ ^ 2Sm^l |c2,m||J2,m|c0s(?72,m - 4'2.,m) ^^^^ 
C0s(?72,o)|c2,o||-F2,o| 

where f]2,m, 4'2,m are the phases of the true CMB and 
the foreground with m = 1 and m = 2 components 
of the quadrupole. Thus, after substitution of k from 
Eq.llSl) into F,q.^, we get 



|C2,o||J2,o| C0s(?72,o) s 
K — [L £), 



|C2,0|COS(C2,0) = |C2,0|COS(772,0) 



1 + Ai' 



(16) 
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where 

M=(|J^2,l|'+l^2,2n/|i^2,0p<l. 

From Eq. p6)) one can get that if e ^ fi, then 
6,0 = '72,0 and |c|'_^ol = A*|c2,o| < |c2,o|- Thus, the 
ILC \c2^o\ mode is smaller than corresponding |c2,o| 
by the factor ^ /i. However, if e 3> /i, we have 

|C2,0|COS(6,0) 

oEm=l |c2,m||i^2,m|cOs(?72 

and now the phase cos(6_o) is determined by the sign 
of the nominator of the right part in Eq. pT|) . This 
result tells us that for some of realizations of the ran- 
dom CMB, the sign of the input signal is no matter. 
The reconstructed ILC phase can be the same 
(^2,o)j or opposite 7^2.0 ± tt to the phase of the true 
CMB. 

To show that this effect takes place in numeri- 
cal simulations of the ILC method, we consider the 
10000 realization input and output maps for the C2,o 
quadrupole component, presented in [3S]. As was 
pointed out in [4] , the difference between the WMAP 
and Eriksen et al. approach [36] for derived CMB sig- 
nals is caused by different choice of disjoint regions 
in the Galactic plane area, rather than the different 
bias. Using an estimator 

Q in out r-r-io r>r~ic: C^Ui 

^ — ^2,0 ' ^2,0 ~ COS ■ COb52,0 i 

where s^q""* = +1 or —1 for the positive and neg- 

j> • r in. out , i- 1 

atlve sign of C2 q components, correspondmgly, we 
have found that for 2148 realizations S — —1. More- 
over, since for the foregrounds ^2 = ^^ir all K-W 
the WMAP bands, practically 43% of the realizations 
having s™q = —I after using the ILC method change 
the sign to s^q = 1 (the effect of flip). We extend 
our analysis to the octupole component of the Erik- 
sen et al. [35] ensemble of input and output signals 
and have found that this effect still takes place, but 
the number of events is slightly smaller than one for 
the quadrupole component. 

In Table 2 we present the number of events for 
S" = — 1, when input and output components of the 
quadrupole and octupole have the different sign. 

One can see from this Table 2, that for the (2,0) 
and (3,1) components the number of events with 
S = —1 reaches the maximum (2148 and 1852 events, 
correspondingly). To show that the flipping of the 
sign of the (2, 0) component can be responsible for 
increasing of the CMB-foreground cross-correlation, 
we take one of the realizations of input Eriksen et al. 
[36] maps, namely "in-00008", and calculated the 
K™{£ = 2) coefficient for this map and the Haslam 



et al. [35] synchrotron map. We have K^'^{i = 
2) -0.222. For this input map, c^^ = -9.944/iA'. 
For the output CMB maps, named as "out-00GG8", 
c§;{,* = 11 AhluK and iC°"*(£ = 2) = -0.522. Then, 
we change the sign of component and recalcu- 
late the coefficient of the cross-correlation once again. 
We have K°'^\t = 2) = -0.2466, which is practically 
the same for the input map. This example clearly 
demonstrates that for input realizations mentioned 
above with negative c^q, this component of the sig- 
nal would be quite likely reconstructed with the op- 
posite sign which leads to increasing a negative level 
of cross-correlation coefficients K°"'*{£ = 2). Roughly 
speaking, the value of the coefficients K depends on 
the flip-effect of the sign for the (2,0) component. 

To show corresponding changes of the images 
of the quadrupole and octupole related to the flip- 
effect, in Fig. [5] we show the ILC (III) quadrupole and 
octupole with the different sign of the (2,0), (2,2) 
and (3,1) components. As one can see the flip-effect 
changes significantly the morphology of the maps. 
This is not surprise since all these £ + m even com- 
ponents are the most powerful components of the 
ILC(III) signal. 

5. CONCLUSION 

We have investigated the cross-correlation of the 
ILC (I) and the ILC (III) with the WMAP MEM fore- 
grounds, and shown that these correlations tightly 
coupled to the bias of the ILC CMB signal. By us- 
ing Monte-Carlo simulations we found the probabil- 
ity distribution function for the coefficient of cross- 
correlation between the true CMB and the foreground 
signals. For debiasing of the ILC CMB, we need to 
know exact value of the coefficient K{£). Since the 
ILC MV method does not provide any additional 
information about the bias, any assumptions about 
the properties of the foregrounds should leave corre- 
sponding residuals on the debiased ILC signal and its 
power spectrum (see for the comparison |H l^ )- 
These uncertainties reflect directly our ignorance of 
the exact value of the realization of the random pro- 
cess K{£), making the reconstruction of the CMB un- 
stable. 

There is another reason for instability of recon- 
struction of the quadrupole and the octupole of the 
CMB related to the CMB-foreground bias. We have 
discovered the flip-effect of the sign for even £ + m 
modes ((2,0) and (2,2) components of the quadrupole 
and the (3,1) component of the octupole). The corre- 
sponding probability of this effect is about 20% (see 
Table 1). This effect could have significant impact on 
the debiasing technique. 
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Figure 5: The illustration of the flip-effect for the ILC(III) quadrupole and octupole. Top is the ILC(III) 
quadrupole, the second from top is the ILC(III) quadrupole with flip-effect of the (2,0) component, the third 
from the top is the quadrupole with flip- effect of the (2,2) component. The second from the bottom and the 
bottom panels show the ILC(III) octupole and the flip-effect for the (3,1) octupole component. 
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